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PRESENTATION OF COMPUTER CODE SPIRALI FOR INCOMPRESSIBLE, TURBULENT, 
PLANE AND SPIRAL GROOVED CYLINDRICAL AND FACE SEALS 


Jed A. Walowit 
Jed A. Walowit, Inc. 
Clifton Park, New York 


OVERALL CAPABILITIES OF PROGRAM 

• Computes rotordynamic coefficients, flow and power loss for cylindrical and 
face seals 

• Treats turbulent and laminar, Couette and Poiseuille dominated flows 

• Fluid inertia effects included 

• Rotordynamic coefficients in 3 (face) or 4 (cylindrical) degrees of freedom 

• Includes effects of spiral grooves 

• User definable transverse film geometry including circular steps and grooves 

• Independent user definable friction factor models for rotor and stator 
User definable loss coefficients for sudden expansions and contractions 
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Coordinate system for seal analysis. 
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ASSUMPTIONS 


• Incompressible and isothermal flow 

. Film thickness small in comparison with other geometric parameters 

• Bulk flow turbulence model 

. Loss coefficients used to treat inertia effects at film discontinuities 

. Axisymmetric primary flow with small perturbations theory used for transient 
and circumferential effects 

• Narrow groove theory used for spiral grooves with inertia treated globally 

• No cavitation 
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BULK FLOW EQUATIONS FOR TURBULENT FLOW 


integrated momentum 


ph 


dU 


du U 5u UVl, 


, +v— + — — + . 

I, at as r ae 


”r ae +(?b_?a) ' T 


ph 



v— u jv 

as r ae 



-h-^P 

as 



integrated continuity 


lA (rvh ) + lJ_ (uh)+ i!l =0 
r as r ae at 



t 

Velocities and forces on a differential element in the 6 direction. 
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SHEAR STRESS AND FRICTION FACTOR RELATIONSHIPS 


^P1U-Q«|J 2HP| “ : — j( G -°.) “ ^ R aMRa)(U-U a ) 




R a = 2h | u - u a | p/ix , Rb 55 2h|u -u b |p/n . 
f a (Ra) = "oFC . W • HofO • Hirs - Blasius 

fa.b = 0.001375 1 + (lO 4 ^ + ^-j 3 - Mood y 



((R,h,v) , Ah < 0 (contraction) 

( i ^ — 1 , Ah 2 0 (expansion) 

1 h-AhJ 
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. _ 1 

Schematic of spiral groove parameters, global and local pressures. 


Print Inputs and check 
for iileoaJ values. Any 
errors? 


| Generate grid. 



y 

Com pul new value of 
C. Are C and axial 
stiffness positive? 


n 


Has Iteration limit i 
been reached? 



i 

i n 




y 

— Set error code 

- 


Compute film and trans- 
verse velocity scale fac- 
tor. Any errors? 

tn ~ 


Compute dimensionless 
variables and parame- 
ters. Any errors? 
fn ~ 


| Set flow direction. | 

Solve first order flow [ 
equations from Inlet to 
exit Any errors? 

. ES 

Is exit pressure con- 1 [1 

dftion satisfied? | " 

R 

Compute load, flow 
and power loss. 

Solve second order flow | 
equations at 0 frequency I 
and at running freqency. 
Any errors 7 

' I" 

Compute stiffness, 
damping and mass 
matrices. 


-I Home In on axial load for 
face seals? 


Print available output], 
and exit. I 


] Set new Inlet velocity. I 


Has Iteration limit 
been reached? 


Flow diagram for overall logic used in computations 
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c CASE 1 ) Cylindrical seal with grooves, laminar, no press, grad. 


&INPUTS 

TITLE 


ELFR 

ALPI 

BETI 

DELT 

ZET 


’Cylindrical seal with grooves, laminar, no press, grad.' 


I FACE 

= 

0 

ISIUN 

*= 

0 




IGROT 

as 

0 

NOI 

* 

2 

I FLOW 

= 

1 

RO 

a 

1 . 0000E+00 

EL 

« 

5.0000E-01 

C 

- 

1.0000E-03 

RPM 

- 

5 . OOOOE+04 

RPMO 

= 

2.5000E+04 

RPMD 


O.OOOOE+OO 

PLEG 

- 

O.OOOOE+OO 

PRIG 


O.OOOOE+OO 

FZD 

55 

O.OOOOE+OO 

vise 

as 

3.0000E-08 

DENS 

= 

O.OOOOE+OO 




EMA 

- 

-2.5000E-01 

ENA 

- 

7.9100E-02 




EMB 

= 

-2.5000E-01 

ENB 

a 

7.9100E-02 




HTAP 

a 

O.OOOOE+OO 

HBRL 

= 

O.OOOOE+OO 




TOLH 

a 

1.0000E-04 

TOLV 


1.0000E-05 

DUT 

- 

1.0000E-06 

I HOME 

m 

0 

NITH 


10 

NITV 

* 

30 

NREG 

- 

2 

NRSUB 

= 

50 50 





5.0000E-01 

5.0000E-01 

2.5000E+01 

2.0000E-03 

O.OOOOE+OO 


5.0000E-01 

0.0000E+00 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 


CYLINDRICAL SEAL, INERTIA NEGLECTED 

LENGTH. DIAMETER . CLEARANCE = 5.0000E-01. 2.0000E+00. 1.0000E-03 (IN) 
ROTOR. SWIRL AND DIST. SPEEDS - 5.0000E+04. 2.5000E+04. O.OOOOE+OO (RPM) 
PRESSURE AT START. END AXIAL BOUNDARIES - O.OOOOE+OO. O.OOOOE+OO (PSI) 
VISCOSITY- 3.0000E-08 (PSI-SEC), DENSITY = O.OOOOE+OO (LB-SEC/IN4) 
ERROR CODE = 0. ITERATIONS IN PRIMARY FLOW = 2 
FLOW = 1 . 1506E+00 (IN**3/SEC) 

TORQUE - 4.3909E-01 (IN-LB), FILM POWER LOSS = 3.4835E-01 (HP) 

AXIAL REYNOLDS NUMBER - O.OOOOE+OO 

CIRC. REYNOLDS NUMBERS FOR ROTOR AT SEAL ENDS - O.OOOOE+OO. O.OOOOE+OO 


DYNAMIC COEFFICIENTS ( FORCE UNIT / DISP. UNIT ) 


DISP. 

Kx 

Ky 

Kphi 

Kpsl 

Bx 

By 

Bphi 

Bpsl 

Ax 

Ay 

Aphi 

Apsi 


x (IN) 
2.7021E+04 
-1.3511E+04 
2.3387E+02 
6 . 5242E+01 
5 . 1297E+00 
-7.8477E-11 
-2.0242E-02 
1.3973E-01 
4.7708E-19 
-4.0552E-18 
1.1927E-19 
-5.4045E-20 


y (IN) 
1.3511E+04 
2.7021E+04 
-6.5242E+01 
2 . 3387E+02 
7.8477E-11 
5.1297E+00 
-1.3973E-01 
-2.0242E-02 
4.0552E-18 
4.7708E-19 
5.4045E-20 
1.1927E-19 


phi (RAD) 
-6.3081E+02 
1 . 2959E+03 
1 . 6699E+02 
-8.6419E+01 
-2.0242E-02 
-1.3973E-01 
3.0437E-02 
-3.8319E-13 
7.4544E-20 
1.7891E-19 
-1.8636E-20 
-3.3545E-20 


psi (RAD) 

- 1 . 2959E+03 
-6.3081E+02 
8.6419E+01 
1.6699E+02 
1.3973E-01 
-2.0242E-02 
3.8319E-13 
3.0437E-02 
-1.7891E-19 
7.4544E-20 
3.3545E-20 
-1.8636E-20 


FORCE UNIT 

LB 

LB 

IN -LB 
IN-LB 
LB -SEC 
LB -SEC 
IN -LB -SEC 
IN -LB -SEC 
LB-SEC2 
LB-SEC2 
IN-LB-SEC2 
IN-LB-SEC2 
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Comparison with results published by D. W. Childs (1983) 




u in = r 0"/2, L/D-1 

u ln =0 - L/D=.2 

i u l n '0, L/D = 1 | 


SPIRALI 

Childs 

SPIRALI 

Childs 

SPIRALI 

Childs 

SPIRALI 

Childs 

Q (cm 3 /$) 

4006. 

4019. 

1771. 

1779. 



1767. 

1779. 

(MN/m) 

18.90 

18.65 

10.79 

9.756 

18.58 

18.52 

13.25 

12.48 

EE 

4.127 

4.213 

91.78 

94.05 

-.3027 

-.3000 



B,;* (KN-s/m) 

21.89 

22.35 

487.2 

500.6 

21.89 

22.47 

489.5 



1.140 

1.206 

102.9 

107.5 

.8518 

.8932 

89.26 

93.39 | 

HS 

3.020 

3.200 

272.6 

285.3 

3.003 

3.189 

272.1 

285.3 | 


DEFINITIONS OF COEFFICIENTS 


Overall Seal Discharge Coefficient 

C d 



Radial Force Coefficient 

-f r = K + co> - Mo 2 = K ef - M ef o > 2 

Tangential Force Coefficient 

-f fl = Co - k = C ef w 
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TITLE 

I FACE 

IGROT 

RO 

RPM 

PLEG 

vise 

EMA 

EMB 

HTAP 

TOLH 

I HOME 

NREG 

NRSUB 

ELFR 

ALPI 

BETI 

DELT 

ZET 


- 'Chi Ids, No! an & Kilgore Stator 1, 25* swirl, 1000 RPM' 

= 0 ISIUN = 1 

= 0 NOI = 0 I FLOW = 0 

= 5 0800E-02 EL = 5.0800E-02 C - 3.5600E-04 

- 1 0000E+03 RPMO = 2.5000E+02 RPMD = O.OOOOE+OO 

= 0.2500E+06 PRIG = O.OOOOE+OO FZD - O.OOOOE+OO 

= 1.5400E-04 DENS ° 1.5700E+03 

= -2.5000E-01 ENA = 7.9100E-02 

= -2.5000E-01 ENB - 7.9100E-02 

- O.OOOOE+OO HBRL = O.OOOOE+OO 

= 1.0000E-04 TOLV = 1.0000E-05 DDT = 1.0000E-06 

= 0 NITO “ 10 NITV = 30 

= 19 

- 5555555555 
555555555 

= .1005 .047 .047 .047 .047 .047 .047 .047 .047 .047 
.047 .047 .047 .047 .047 .047 .047 .047 .1005 

- 1. 1. 1. 1. 1. 1. 1. I- I- I- 

1. 1. 1. 1. 1. 1. 1- 1- 1. 

= 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 

0. 0. 0. 0. 0. 0. 0. 0. 0. 

= 0. 3.8E-4 0. 3.8E-4 0. 3.8E-4 0. 3.8E-4 0. 3.8E-4 0. 3.8E-4 
0. 3.8E-4 0. 3.8E-4 0. 3.8E-4 0. 

- 0 . 0 . 0 . 0 . 0 . 0 . 0 . 0 . 0 . 0 . 

0. 0. 0. 0. 0. 0. 0. 0. 0. 



A P (bar) 


Dimensionless axial flow rates 


16 





mm r 
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